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Abstract 

Starting from general considerations about the different concepts 
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1 Keeping up with experiments 



Besides existing machines and a whole bunch of data sets, further lepton pair 
coUiders such as the Next Linear CoUider (NLC) and the Future Muon Col- 
Uder (FMC) are designed for the future. They will be able to determine the 
properties of heavy particles to an extend never reached by present colliders. 
Especially they will be able to measure the properties of the top quark which 
was first measured at the Tevatron with a mass of rrit = 174.3 ± 5 GeV (see 
e.g. |1[ and references therein). Although the top quark will be studied at 
the LHC and the Tevatron (RUN-II) as well with an expected accuracy for 
the mass of 2 — 3 GeV, the measurements at NLC are expected to have an 
accuracy of 0.1% (100 - 200 MeV) §. 

Theorists are therefore asked to improve their predictions, too. Of spe- 
cial interest is the consideration of the cross section near the pair production 
threshold because the resonance observed there is constituted solely by a 
colour singlet state. A lot has been done in the past and just recently to cal- 
culate and match corresponding diagrams up to the next-to-next-to-leading 
order in order to determine this cross section (for an overview on work that 
has been done see e.g. Ref. 0). But the question remains what we actually 
want to determine. What can we say about the mass of a particle which 
never will show up in nature in isolated fashion because of its confinement? 
Does the term of the mass of the top quark really makes sense? 

The everydays life of theorists show that this actually makes no sense. 
There are different concepts on the market - each of them with different 
advantages and disadvantages, regions of validy and failure. They all show 
that there cannot be a unique definition of the quark mass. Instead the task 
of the theory can only be to quantify and limit the theoretical uncertain- 
ties to an amount comparable to the experimental measurements. One of 
these theory-born uncertainties are special kinds of IR singularities known as 
renormalons which limit the accuracy of a prediction to the order of Aqcd- 
Just recently this problem has been overcome (a synapsis of the progress is 
given in Ref. and this talk will deal with the different concepts which 
were developed in this context, closing with an own suggestion presented by 
the authors 0. 
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2 MS mass and on-shell or pole mass 



A mass quite often used in perturbation theory is the MS mass. It is related 
to the bare mass of the underlying QCD by the calculation of self energy 
diagrams of the quark. The term MS mass indicates that for this calculation 
we use the MS subtraction scheme which itself is based on the dimensional 
regularization method. The calculation of self energy diagrams which is done 
actually up to the third order QCD perturbation theory 0, ^ leads to a 
couple of master diagrams which are calculated at center-of-mass energies in 
the deep Euklidean domain (g^ < 0, |g| large) where the perturbation series 
converge sufficiently good. This MS mass of course does not fit the range 
around the threshold and spoils the non-relativistic expansion if applied. But 
the MS mass can be continued to this region. 

The continuation is done by the Borel transformation. This transform 
can be understood as the extrapolation from the deep Euklidean domain to 
the tt threshold domain where a simultaneous limit of the negative center-of- 
mass energy and the order of the used derivative is taken. The ratio of these 
two quantities remains fixed and constitutes the Borel parameter. As shown 
in Ref. [0], this Borel transformation can be made explicit for the case of the 
"large /3o limit" where the prescription of Brodsky, Lepage and Mackenzie p| 
is used to sum up all fermion loop insertions (being proportional to the zeroth 
order QCD beta function coefficient Pq = 11(7^/3 — 2 A^'^^^/S) in order to obtain 
an estimate (for a relation of this estimate to the MS mass see Ref. 0). 

If we use the prescription for the static quark- ant iquark potential V{r), 
we obtain 

oo 

Vir) = E VniPoasT^'. (1) 

n=0 

The Borel transformation of this series expression leads to 

V{r-u) = f:^u- (2) 

The series for the Borel transformed quantity converges even if the pertur- 
bative series for V{r) itself is divergent. The inverse Borel transformation 

POO _ 

V{r) = / e"/(^o°^V(r;M)rfM (3) 
Jo 

provides then a good definition. However, there might occur singularities on 
the integration contour. As mentioned before, for the "large /3o" estimate 
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(also known as "naive non-abelianization" ) we obtain an explicit expression 



~ V{l/2 + u)V{l/2-u) 
^(^'"^ f{2^) 

which has simple poles for u = k + 1/2, k = 0,1, 2,.... These singularities 
are a special case of IR renormalons. 

The same happens then for the continuation of the MS mass to the on- 
shell mass in the threshold region. This on-shell mass is appropriate to the 
treatment near the threshold because the quark is then assumed to be nearly 
on-shell. The on-shell mass m (or pole mass as it is called more often) can be 
determined by the demand that the inverse Fermi propagator should vanish 
at this point. 

But the continuation of the MS mass via Borel transformation to the pole 
mass leads again to IR renormalon ambiguities. 

So there is no way out? Wait, there is one! Both the pole mass and the 
static potential show an IR renormalon ambiguity. Maybe they cancel in a 
specific combination of these two quantities? Indeed it emerges that in the 
combination 2m + V{r) the IR sensitivity vanishes. Only for this reason it is 
then possible to construct a quantity which has no such ambiguities. In fact 
there are several of them, called threshold masses which we will deal with in 
the following. 



3 Threshold mass definitions 

The first definition we want to mention here is the definition given by Hoang 



and Teubner [|10[. We mention it first because it is the most sophisticated 



definition. Starting with the non-relativistic QCD (NRQCD) [|rT|, an effec- 
tive theory designed to handle the non-relativistic quark-antiquark system 
which separates the low-momentum scales mv and mv'^ from the high mo- 
mentum scale m {v is the velocity of the quark), they stress the relevance 
of the different scales v and ag and their interrelation. Especially in the 
non-relativistic regime even the space and time components can be scaled 
differently. So there are "soft" (g ~ m{v;v)), "potential" (g ~ m{v'^;v)), 
and "ultrasoft" regions (g ~ m(t>^;t>^)) where the last one is populated only 
by gluons, not by quarks. If we separate the dynamical degrees of freedom 
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also from the "soft" heavy quarks and gluons and the "potential" gluons, 
supplemented by an expansion in momentum components to the order mf ^, 
this leads to the construction of the "potential NRQCD" (PNRQCD) [g 
with the Langrangian 

>CPNRQCD = -Cnrqcd + / {i^HW)V {r){x\){^)d?r. (6) 

Here the heavy quark- antiquark potential V{r) occurs, and x ^^^e the quark 
and antiquark spinors, respectively. Finally the authors of Ref. come to 
the conclusion that an adequate threshold mass definition is given by one 
half of the perturbative mass of a fictions toponium l^S*! ground state. This 
threshold mass definition is referred to in the literature as IS" mass. 

A second threshold mass definition is given by Bigi et al. [1T3 and is known 
by the name "low scale" (LS) mass. 

Finally, Beneke gave a definition which then led to the generalization sug- 



gested by us. The static potential subtracted (PS) mass [0 is the minimal 
version of a threshold mass. As the name tells us, the definition is given by 
combining IR singular parts taken from the static quark-antiquark potential 
with the pole mass in order to cancel out the IR sensitivity. Starting point 
is the observation that the IR sensitive parts of the potential stem from long 
distance contributions, i.e. small values of the three-momentum. If we there- 
fore restrict the three-momentum range of the potential as expressed by the 
Fourier transformation 

V{\r\) = J (03^^'''^(l^1) (7) 

by |g| < Aqcd, this results (for the leading order Coulomb potential) in the 
contribution 

L, « Jo:^^ ^ ~ 7^ ^QCD + (^(Aqcd^ j (»j 

which restricts the accuracy. The way out is to use a potential V{r;^f) 
defined by the Fourier transformation with |g| > /ij and to subsume the 
remaining part to the mass, 

1 f dPq ~ 

mpsifif) =m-6mps{f^f), 6mps{fif) = -- TT^^d'?!) (9) 

Z J\q\<^lf [ZlT)-^ 

(it is legitimate to use e*'''^' ~ 1 here). The scale /i/ introduced by this is 
called factorization scale. It drops out again if the static PS mass is related 
to the MS mass. 



5 




g |5 



P Pi Pn P 

Figure 1: The general structure of the self energy diagram of a quark 



4 A dynamic mass definition 



It has been pointed out already in Ref. [|T4| that the replacement 

^ -i7r6{v ■ q) (10) 



V ■ q 



in the self energy leads to the leading infrared behaviour. We have taken this 
consideration as starting point for our definition and extended it to the more 
general non-static case. 

The general structure of the self energy diagram for an on-shell quark 
with mass m and momentum p {p^ = m^) is shown in Fig. |I|. The expression 
for this diagram is given by 

m=l \^^) 

where the last factor is a non- commutative product with decreasing index 
n. The line momenta fc„ are linear combinations of the gluon loop momenta 
Im, the particular representation is specified by the structure S. The symbol 
{Im} means the set of all these loop momenta, the same symbol is used for 
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the Lorentz and colour indices. In general we have M < N which means that 
line momenta can be correlated. The momenta of the virtual quark states 
are given by p„ = p + /c„. 

The quark is considered to be at rest, p = (m, 0), and it interacts with a 
number of gluons. The subdiagram 5* displayed in Fig. |l| describes the inter- 
action between the gluons. In general the quark lines between the interaction 
points represent virtual quark states. However, if the virtual quark comes 
very close to the mass shell and the total momentum of the cloud of virtual 
gluons becomes soft, this situation gives rise to long-distance nonperturba- 
tive QCD interactions. The described (virtual) contributions results in the 
soft part of the self energy, Sgoft- So we define 

N M 74/ 

-^^soM = E / n J^AZlii^rn}) {-^9sr--^T^.,.) X 

i=l m=l l^^J 
i+l ^ 

^ n / i-igsl'^-TaJ + m) [-iTt6{p^i - m^)) x 

n=N fri ^ 

X (-z<7.7"'T,J n -r^{-i9sl''"Ta:). (12) 
ntt-i V>n-m 

One can derive this expression from Eq. ([TT|) by using the identity 

= _ _^ p (13) 

p^ — + It yp-^ — m'^ J 

and the fact that the principal value integral does not give any infrared sen- 
sitive contribution. The delta function can be used to remove the integration 
over the zero component of fcj. In order to parameterize the softness of the 
gluon cloud we impose a cutoff on the spatial component, \ki\ < /ij, and 
indicate this by the label fif for the factorization scale written at the upper 
limit of the three-dimensional integral. So we can rewrite Eq. (p!^ ) as 

Ssoft(i^,/x/) = -lip (14) 

where 

„Af-l Ml 

Vik,p) = - li 7^^{t}({^-}) (-^^?^7"'^"^^<^^^+0 X (15) 

m=l l^^J 

n=N Pn-m 2p'> ^t^Li Pn-m 
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and finally define 



rrtpg = m — 5mpg with 5mpg 



^soft I 



(16) 



In the following we deal with the different realizations of this compact expres- 
sion. As we will see explicitly, the function V{k,p) occurring as integrand can 
be seen as quark- ant iquark potential where we have summed over the spin 
of the tensor product of a final state with an initial state. Because the static 
quark- ant iquark potential is used in a similar way in Ref. 
result of Ref. WM in the static limit. 



14 , we recover the 



5 The one-loop contribution 

The leading order contribution to the self energy of the quark is given by 

where Feynman gauge is used for the gluon. The soft part of it is given by 
S,oft(#,/i/) = --y (^^(^'P)' V{k,p) = V+{k,p) + V.{k,p) (18) 
where 



V^{Kp) = 9lC,Wm^ + .^T2m)_ ^^^^ 

2m\/iv? + K^{\/im? + — m) 

and K = \k\. The procedure of taking the soft part by setting a cut is shown 
in Fig. ^(a-b). The two potential parts are known as the scattering potential 
Vj^{k,p) and the annihilation potential V_{k,p) and correspond to the two 
zeros k± = ±v^/?~+^r? — m of the delta function. Integrating the potential 
up to the factorization scale fJ^f, we obtain 



E.o,,(^,) = 3in + + 1 - + 1 ^ (20) 



The expansion of this expression in small values of fJ^f/m results in 

fi_^,o(^]|. (21, 



<X < ^ <X <x 

(a) (b) (c) (d) 

Figure 2: Leading order contribution to the quark self energy (a) and to the 
quark- ant iquark potential (b). The cross indicates the point where we cut 
the quark line by imposing an on-shell condition to the virtual quark state. 
In Coulomb gauge the gluon propagator can be decomposed in a Coulomb 
propagator (c) and a transverse propagator (d). 



The first term reproduces the result given in Ref. to leading order in 
as while the second term is the recoil correction to the static limit in this 
order of perturbation theory. This second term is related to the Breit-Fermi 
potential but does not coincide with it. 



6 Two-loop contributions 

To take a step beyond the leading order perturbation theory, we consider 
two-loop diagrams for the heavy quark self energy as shown in Fig. ^. In 
cutting the quark line in all possible ways we obtain a lot of diagrams from 
the ones shown in Fig. |^. However, we find that the final contribution of 
the two abelian diagrams in Fig. ^(a-b) to the soft part of the self energy is 
suppressed by fi'j/m'^. Therefore it turns out that the only abelian diagrams 
which can give a non-suppressed contribution to the soft part of the quark 
self energy are the diagrams containing the vacuum polarization of the gluon 
as shown in Fig. ^(d-f). The simple calculation of these diagrams within the 
MS scheme, accounting only for light fermion loops, gluon loop (and ghost 
loop if Feynman gauge is used) results after renormalization in 



(22) 
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with the same coefficient ai as in Ref. fl^] (see below). This result has 
been anticipated because the expression in the curly brackets of the inte- 
grand reproduces the next-to-leading order correction to the QCD Coulomb 
potential. 

For the only non-abelian diagram in Fig. ^(c) we obtain 

soft ~ /^/- V^"^/ 

This result has been anticipated, too, to be minus one half of the non-abelian 
correction to the QCD Coulomb potential, which is known in the literature 
(see for example Refs. [^, 0), 



7 Our final result 

Summarizing all contribution up to NNLO accuracy, we obtain 

10 



(a result to a not yet specified order in Hf/m) where m is the pole mass, ix 
is the renormalization scale, /ij is the factorization scale, and 



Co — 1, Cq — 0, Cg — 



2' 

2 



Ci = ai-2/3oln(^^j , C( = C^y, 

C2 = a2-2(2ai/3o + A) (^ln(^^^ -1^ + 

+4,„^(ln^(/;i)-21n 2). (26) 

The constants ai and 02 are given by [0, [19 



ai = -Ca--TfNf, 

(AM?, . vr^ 22 \ ^2 56 



V 162 4 3 ^7 ^ V 81 

+ f ^T^iVp) ' - f ^ - 16(3) CfTfNf, (27) 



9 / V 3 

the coefficients of the beta function are known as 

11 4 34 2(1 

= -Ca - -TfNf, (3, = -Cl - -CaTfNf - ACfTfNf (28) 

where Cf = 4/3, Ca = 3 and Tp = 1/2 are colour factors and Nf = 5 is 
the number of light flavours. The coefficients Ci and C2 in Eq. (^5]) have 
been derived in Ref. [|T1| by using known corrections to the QCD potential. 
In Ref. we have obtained the coefficients Cg, Cg, and C[. Note that our 
result can be represented in a condensed form as 

m^{^^f) -^=-\r {ycm + Vnm) + V^Am) (29) 

where the first term Vc is the static Coulomb potential, Vr is the relativistic 
correction (which is related to Breit-Fermi potential but does not coincide 
with it), and Vata is the non-abelian correction. Note that we did not include 
electroweak corrections in our calculations. 
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m{m) 




™NLO 


^NNLO 
// ipg 

™NNLO 




™NLO 
'"pole 


™NNLO 
"'pole 


160 


166.51 


167.69 


167.97 
167.95 


167.44 


169.05 


169.56 


165 


171.72 


172.93 


173.22 
173.20 


172.64 


174.28 


174.80 


170 


176.92 


178.17 


178.47 
178.45 


177.84 


179.52 


180.05 



Table 1: Top quark mass relations for the MS, PS, PS, and the pole mass at 
LO, NLO, and NNLO in GeV. We fix the MS mass to be m(m) = 160 GeV, 
165 GeV, and 170 GeV and find the pole mass at LO, NLO, and NNLO from 
the three-loop relation in Refs. ||^, The PS and PS masses are derived 
from the pole mass by using Eq. ( P^D (without or with the 1/m contributions, 
respectively). We use the QCD coupling constant agimz) = 0.119 pO| , 



m[m] 



and /i/ = 20 GeV. 



Using the relation given by Eq. (^3]) between and m as well as the 
relation between m and mijn) given in Refs. H, we fix the MS mass to 
take the values mini) = 160 GeV, 165 GeV, and 170 GeV and determine the 
pole mass at LO, NLO, and NNLO. This pole mass is then used as input 
parameter m in Eq. (^Sf ) to determine the PS and PS masses at LO, NLO, 
and NNLO. The results of these calculations are collected in Table |I|. Taking 
the PS mass instead of the PS mass, we observe an improvement of the 
convergence. The differences for the mass values in going from LO to NLO 
to NNLO for e.g. m(m) = 165 GeV read 7.64 GeV, 1.64 GeV, and 0.52 GeV 
for the pole mass, 6.72 GeV, 1.21 GeV, and 0.29 GeV for the PS mass and 
finally 6.72 GeV, 1.21 GeV, and 0.27 GeV for the PS mass. So the tendency 
is going in the proper direction. However, we stress that the goal, namely to 
avoid the infrared renormalon ambiguities, is already reached by the static 
PS mass as well as by the other threshold masses which have been developed. 
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8 What to do with our result then? 



Having mass values for the top quark at hand, we can consider the cross 
section of the process e^e~ — > tt in the near threshold region where the 
velocity v of the top quark is small. It is well known that the conventional 
perturbative expansion does not work in the non-relativistic region because 
of the presence of the Coulomb singularities at small velocities v 0. The 
terms proportional to (as/f)" appear due to the instantaneous Coulomb 
interaction between the top and the antitop quark. The standard technique 
for resumming the terms (ag/v)"' consists in using the Schrodinger equation 
for the quark-antiquark potential 



{H-E- iT)G{r, r'\E + iV) = 6^ 



(30) 



where H is the non-relativistic Hamiltonian of the heavy quark-antiquark 
system, to find the Green function. The Green function is then related to 
the relative cross section by the optical theorem, 

ti) 



R 



(T(e+e~ - 
2 .r 247r 



C(ro)Im 



1 - 



3m^ 



(31) 



ro-O 



where a{e~^e^ A^^A* ) = 47ra^/3s, eg is the electric charge of the top 
quark, Nc is the number of colours, y/s = 2m + E is the total center-of-mass 
energy of the quark-antiquark system, m is the top quark pole mass and F 
is the top quark width. The unknown short distance coefficient C(ro) can 
be fixed by using a direct QCD calculation of the vector vertex at NNLO in 



the so-called intermediate region ||22| , P^] and by using the direct matching 



procedure suggested in Ref. p^ . 

If we take the PS mass as input parameter instead of the pole mass, the 
energy value E has to be changed to ^i — 2mps(/i/) instead of y/s — 2m. The 
same holds for the PS mass. It is shown that the Schrodinger 

equation can be reduced to the equation 

{Hi- Ei)Gi{r,r'\Ei) =6\r -r') (32) 

with the energy Ei = E + E^/Am, E = E + iT, and with the Hamiltonian 

3E,_ . /2 C/ 



Hi = — + Vcir) + ^Voir) 
m 2m 



3 Cf. 



2m 
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as(At)C'F 

«.(/^)CF^2/5oln(/iV) + ai) + (33) 



Ait 



/?2 |^41n2(/iV) + y) + 2(/5i + 2/3oai) ln(/iV) + | 



where /i' = /xe'*'^, /i is the renormahzation scale, and je is Euler's constant. 
The final expression for the cross section is then given by 



R = ^C(ro)Im 



1-1^) Gi(ro,ro|i?i) 



(34) 



6m ^ 

with Gi(ro, rol-Ei) being the solution of Eq. (152) at r = r' = tq. For the 



numerical solution we use the program derived in Ref. |2^. Note that we 
do not take into account an initial photon radiation which would change the 
shape of the cross section. This can be easily included in the Monte Carlo 
simulation. 

Taking the pole mass as input parameter, the top quark cross section 
at LO, NLO, and NNLO is shown in Fig. ^ as a function of the center-of- 
mass energy. For the top quark pole mass we use rrit = 175.05 GeV, for 
the top quark width Fj = 1.43 GeV, and for the QCD coupling constant 
as{mz) = 0.119 0. Different values fx = 15 GeV, 30 GeV, and 60 GeV for 



the renormahzation scale are selected because they roughly correspond to 
the typical spatial momenta for the top quark. We see that the NNLO curve 
modifies the line shape by the amount of 20 — 30% which is as large as the 
NLO correction. It also shifts the position of the 15* peak by approximately 
600 MeV. These large shifts of the peak position are expected because of the 
renormalon ambiguity. 

If we do the same for the PS mass value mps(/i/ = 20 GeV) = 173.30 GeV 
which corresponds to the pole mass m = 175.05 GeV according to the "large 
Pq" accuracy estimate, we obtain a picture like in Fig. ^ where we observe an 
improvement in the stability of the position of the first peak in comparison to 
the previous analysis as we go from LO to NLO to NNLO. Actually, for the 
three values fi = 15 GeV, 30 GeV, and 60 GeV we obtain the maxima of the 
IS peak for NNLO at w = 347.32 GeV, 347.41 GeV, and 347.48 GeV while 
the maximal values are given by -Rmax = 1.379, 1.184, and 1.088, respectively. 
This demonstrated that a large variation in the renormahzation scale fi gives 
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Figure 4: The scheme using the pole mass: shown is the relative cross section 
R{e^e~ ^ tt) as a function of the center-of-mass energy in GeV for the LO 
(dashed-dotted lines), NLO (dashed lines), and NNLO (solid lines) approx- 
imation. We take the value rrit = 175.05 GeV for the pole mass of the top 
quark, Fj = 1.43 GeV for the top quark width, as{mz) = 0.119 and different 
values /X = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 



rise only to a shift of about 160 MeV for the IS* peak position at NNLO 
while the variation for -Rmax is still large. The same procedure done for the 
PS mass shows no observable difference. The differences shown in Table || 
are too small to make any viewable effect. 



Conclusion and discussion 

We discussed the potential subtracted (PS) mass suggested in Ref. in 



clS db 



definition of the quark mass alternative to the pole mass. In contrast to the 
pole mass, this mass is not sensitive to the non-perturbative QCD effects. 
We have derived recoil corrections to the relation of the pole mass with the 
PS mass. In addition, we demonstrated that, if we use the PS or the PS 
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Figure 5: The scheme using the PS mass: shown is the relative cross section 
R{e^e~ ^ tt) as a function of the center-of-mass energy in GeV for the LO 
(dashed-dotted hues), NLO (dashed hues), and NNLO (sohd hues) approx- 
imation. We take the value mps = 173.30 GeV for the PS mass of the top 
quark, Fj = 1.43 GeV for the top quark width, as{mz) = 0.119 and different 
values /i = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 



mass in the calculations, the perturbative predictions for the cross section 
become much more stable at higher orders of QCD (shifts are below 0.1 GeV). 
This understanding removes one of the obstacles for an accurate top mass 
measurement and one can expect that the top quark mass will be extracted 
from a threshold scan at NLC with an accuracy of about 100 — 200 MeV. 
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